Complexity spectrum of some discrete dynamical systems 



N. Abarenkova 

Centre de Recherches sur les Tres Basses Temperatures, B.P. 166, F-38042 Grenoble, France 
Theoretical Physics Department, Sankt Petersburg State University, Ulyanovskaya 1, 198904 Sankt Petersburg, Russia 

J.-Ch. Angles cTAuriac 
Centre de Recherches sur les Tres Basses Temperatures, B.P. 166, F-38042 Grenoble, France 

S. Boukraa 

Institut d'Aeronautique, Universite de Blida, BP 270, Blida, Algeria 
J.-M. Maillard 

OO ■ LPTHE, Tour 16, ler Stage, 4 Place Jussieu, 75252 Paris Cedex, France 

Q\ ' (February 5, 2008) 

o\ : 

We first study birational mappings generated by the composition of the matrix inversion and 
of a permutation of the entries of 3 x 3 matrices. We introduce a semi-numerical analysis which 
enables to compute the Arnold complexities for all the 9! possible birational transformations. These 
complexities correspond to a spectrum of eighteen algebraic values. We then drastically generalize 
£NJ ■ these results, replacing permutations of the entries by homogeneous polynomial transformations of 

04 ' the entries possibly depending on many parameters. Again it is shown that the associated birational, 

or even rational, transformations yield algebraic values for their complexities. 
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C*~) . Key words : Arnold complexity, discrete dynamical systems, rational mappings. 

o 
o 

OO 

On 

*3 



O 
c3 



I. INTRODUCTION AND RECALLS 



Birational transformations have been seen to be a powerful tool to analyze the symmetries of the parameter space 
of lattice models of statistical mechanics [^jQ] and to seek for some possible Yang-Baxter integrability Beyond 
the lattice statistical mechanics framework, birational transformations are worthy to be studied per se, as discrete 
dynamical systems. Discrete dynamical systems have been intensively studied (see for example [^[l)- Among them 
polynomial examples, like the Henon map 0, have been precious to understand some features of chaos. Beyond, 
rational mappings are of special interest since they allow some analytical calculations. Furthermore, the rational 
transformations also allow numerical calculations which can be performed with any wanted precision : the existence 
of singularities in the rational transformations one iterates, and their possible "proliferation" is not in fact a numerical 
. ,-h , obstruction. We will first consider mappings generated by the composition of the matrix inverse and some arbitrary, 
but fixed, permutation of the entries of q x q matrices. The results, displayed in this paper, are given for q = 3, but 
are actually valid, mutatis mutandis for arbitrary q values. 

_ _ i 

A. Recalling a previous 3x3 analysis 

Integrability of a mapping amounts to saying that all the orbits of the iteration correspond to elliptic, or rational, 
algebraic curves ||,|9| . From the point of view of the growth^ of the complexity of the successive iterations |^J^] , such 
integrability in curves always yields a polynomial growth of the calculations |0,0], instead of the exponential growth 
one generically expects. Conversely polynomial growth is not restricted to integrability in curves but may correspond 
to orbits "densifying" Abelian varieties 



1 When one iterates a rational transformations the "size" of the successive rational expressions, corresponding to the iV-th 
iterate, grows, in general, exponentially. In particular the degree of these successive rational expressions has, generically, an 
exponential growth jj],[to]]. Growth of the calculations related with factorizations were also introduced by Veselov for some 
particular Cremona transformations Jl l| — |l 3|] . 
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A first exhaustive analysis of all the 9! birational transformations generated by the composition of the matrix 
inversion and of a (fixed) permutation of the entries of 3 x 3 matrices has already been performed concentrating on 
the extraction of integrable mappings |Q. This analysis was exhaustive, but restricted to particular integrability 
criteria]^. Even from this "integrability-digger" point of view some integrable mappings are missing (for example the 
so-called "Class III" mappings of fla], as well as some polynomial growth situations). In the first part of this paper 
we will revisit these 9! — 362880 birational mappings without any a priori integrability criterion and with the help of 
a new equivalence relation among permutations (symmetry) . This analysis exactly yields all the polynomial growth 
situations, and, far beyond, classifies the exponential growth situations. The classification relies on the value of the 
Arnold complexity ]l6[ of the mapping. This complexity can be obtained jlTj from generating functions associated 
with factorization schemes HI detailed below. 



B. Factorization scheme and generating functions 

We use the same notations as in |^|^,[l5|] , that is, we introduce the following two transformations, the usual matrix 
inverse I and the homogeneous matrix inverse / : 

T: M — ► M -1 , and: I: M — > det(M ) • M ~ 1 (1) 

The homogeneous inverse I is a homogeneous polynomial transformation, which associates, with each entry of Mo, 
its corresponding cofactor. Transformation t is any (fixed) permutation of the entries of the 3x3 matrix. We also 
introduce the (generically infinite order) transformations : 

K = t ■ I and K = t-T (2) 

Transformation K is clearly a birational transformation |^,^|. 

For all the various birational transformations associated with permutations of the entries of 3 x 3 matrices, the 
following factorization relations happen to occur at each iteration step Q : 

/i = det(M a ) , Mi = K (M ) , f 2 = — , M 2 = m , h = , , , M 3 = m Vo , 



ft ' fx ' ft ■ ft fr-f? 



and for arbitrary n 



det(M„) = f n+ x-ft-fti-ft2-fZU-ftU---ft (3) 
K(M n ) = M n+1 - f r r flU- fT^- flU- flU--- ft- 1 (4) 
det(M„) • M n+1 = (fl" +1 • • flU ■ flU ■ flU ■ ■ ■ /f") • K(M n ) (5) 

defining the positive integer exponents rj n , (f> n and p n . The / n 's are homogeneous polynomials of the entries of M . 
These factorizations allow to define, at each iteration step, the successive / n 's one can "factor out", and the "reduced 
matrices" M n 's, such that their entries are homogeneous polynomial expressions of the initial entries, and have no 
further factorization. One finds out, looking at the first thirties iterations, that one recovers the same exponents (r) n , 
4> n , Pn) at each iteration step (up to the last emerging coefficient for f\). We assume that this regularity property 
holds for arbitrary n. This regularity^ property assumption is crucial in our analysis. 

We will denote a n the degree of the determinant of matrix M„, and (3 n the degree of polynomial /„ and a(x), 
(3(x), n(x), (j)(x) and p(x), the generating functions of the degrees a n 's, fin's, and of the exponents ?7n's, p n 's and <fi n 's 
in the factorization schemes : 



a ( x ) = ^a n -x n , (3(x) = -x n , r)(x) = }^ r\ n ■ x n , </>(x) = 0„ ■ x", p(x) = y^p n 

n— n— 1 n— n— 1 n— 



2 Associated with particular recursions on some "determinantal" variables XnS, also introduced here in section (tB). 

3 In fact it is shown in |[^] that other slightly more general factorizations scheme can occur on some if-invariant subvarieties 
(yielding smaller Arnold complexity values). Such slightly more general factorizations scheme will also be detailed below (see 
Appendix B). However for the transformations K associated with permutations of q x q matrices, for a generic initial matrix, 
one gets factorization schemes like (N), M), also depicted in |Q|. 
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where ao = 3 and pi = 3. It is straightforward to show ^ that the existence of the stable factorization scheme (||), 
(||) yields the following simple linear relations between these various "degree generating functions" and "exponents 
generating functions" : 

a(x) + 3 • x ■ n(x) • P(x) = 3 + 2 • x ■ a(x) (6) 
x ■ a(x) — 4>{x) ■ P(x) (7) 
3 + 3 • p(x) ■ /3(x) = (l+x)-a{x) (8) 

When analytically iterating an arbitrary transformation K, the degree of the successive polynomial expressions one 
encounters, grow exponentially : a n or (3 n ~ A™. A measures the grows of the calculations and identifies with the 
notion of Arnold complexity |l6| , |l7j ]. From now on A will be called the complexity. When the degree generating 
functions a(x) or P{x) happen to be rational functions, the complexity A is obviously the inverse of the pole of 
smallest modulus. Recalling the "determinantal" variables ||[){l5| x„'s defined by : 



x n {M ) = det(K n+i (M ))-det(K n (M Q )) (9) 

one finds out that these determinantal variables happen to decompose on a product of the homogeneous polynomials 
/„'s only : 

x n (M Q ) = fZi • IT ■ fT-i ■ fT-2 ■ /"-a ■ fT-A ■ ■ ■ (10) 

which defines some, at first sight, "new" exponents^] u>„'s and consequently a, at first sight, "new" generating function 
W(x) : 

OO 

W(x) = ^2w n -x n (11) 

It is worth noticing that the determinantal variables x„'s induce the homogeneous polynomials / n 's emerging from 
the factorization schemes (|]), ([|) and no other homogeneous polynomials. The variables x„s are well-suited since 
they are invariant under a multiplication of Mq by a constant : Mq — > Cst • Mq. In other words the x ra 's are 
homogeneous expressions of degree zero. Concentrating on the degrees of the left-hand side, and right-hand side, of 
(|l0|), one gets the following "degree equation" : 

= f3 n+ l ■ W + /3 n ■ Wi + • • • + p n - p ■ W p+ i + ■ ■ ■ + Pi ■ w n (12) 

from which one immediately deduces the simple functional equation : 

W(x) ■ P(x) = Pi ■ w = 3 • x (13) 



This result is immediately generalized to q x q matrices. Relation ([L3|) becomes W(x) ■ P(x) = q ■ x. From (|13j) 
one actually sees that W(x) is not a new generating function : it is simply related to the degree generating function 
P(x). The complexity A is associated to the zeroes of W(x). 

From relations @, (0) and (||), one easily gets the "degree generating functions" a(x) and P(x) from two of the 
"exponent generating functions" (for instance <f>{x) and rj(x) or rj{x) and p{x)). As a matter of fact, for most of the 
permutations, the factorization schemes are periodic (r) n = ?7„+at, <p n — 4> n +N and p n — p n +N for some integer N). 
Consequently, the exponent generating functions 4>(x) and rj{x) , or p(x), are rational functions with A^-th root of 
unity poles [g| (see, for instance, the exponent generating function p(x) in ( |i~4] ) or (^6|)). In a second step one deduces, 
from relations (|6|), or (0), rational expressions for the degree generating functions a(x) and P(x). However it will be 
seen below that, for some permutations, the factorization schemes are still regular, but not with periodic exponents 
(see the exponent generating function p(x) in ( |l5| ) or (p7|)) : the exponents rj n , cf> n and p n grow exponentially, but one 
remarks that the associated generating functions rj(x), <fi(x), p{x) are still rational, and thus a(x) and P(x) are also 
rational. The exponent generating functions can be seen as an "encoding" of the degree generating functions, and 
thus of the complexity A. Remark that all these rational expressions involve integer coefficients, yielding algebraic 
values for their poles and for the growth of the calculations : the degrees of the successive rational expressions, namely 



4 The Wn's are relative integers and not natural integers like exponents 77™, 4>n and p n 
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a„'s and /3„'s grow like A™, where A is an algebraic number, and for regular factorization schemes (like (15) or (p.7|)), 
see below) , the exponents rj n and <f> n grow like fx" , where fi is the "scheme complexity" . Note that \i is obviously 
such that [i<\. Exponent fj, is also the inverse of the pole of smallest modulus of the exponent generating functions. 
Complexity A allows all kinds of handy, efficient, and formal, or semi- numerical, calculations. We will present below 
such a semi-numerical method and apply it to all the permutations of entries of 3 x 3 matrices^. 

II. COMPLEXITY SPECTRUM ANALYSIS FOR PERMUTATIONS 
A. A semi-numerical method 

All these considerations allow us to design a semi-numerical method to get the value of the complexity A for the 
iteration of rational transformations. The idea is to iterate, with K, a generic initial matrix with integer entries. 
After one iteration step the entries become rational and we follow the magnitude of the successive numerators and 
denominators. During the first few iteration steps some "accidental" simplifications may occur, but, after this tran- 
sient regime, the integer denominators (for instance) grow like A™, where n is the number of iterations. One can 
systematically improve the method as follows : the initial matrix is chosen in such a way that it avoids, as much 
as possible, any "accidental" additional factorization in comparison with the factorization scheme associated with a 
generic matrix. For instance, in a factorization scheme framework like (^), (|J), one chooses the initial matrix Mo 
with integer entries such that the determinant, and most of its cofactors, are prime numbers as large as possible. 
One may impose further constraints on the initial matrices, for instance, that the first homogeneous polynomials fa 
and /a are also, as large as possible, prime numbers. These conditions down-size the probability that all the entries 
of the reduced matrices M n , or the polynomials /„'s, could be divisible by some accidental additional fx, f 2 or / 3 . 
Such initial matrices, well-suited for the iteration of the homogeneous transformation K, are also well-suited for the 
iteration of the (bi-)rational transformation K. In practice we start with a set of initial matrices and keep only the 
one for which the less factorizations occur (non-generic factorization can only correspond to additional factorizations). 

The computations are done using an infinite precision C-library [ jl8| . We perform as many iterations as possible 
during a given CPU time T. This number of iterations, n, is such that T ~ A™. For A close to 2 and T = 60 
seconds, n is of the order of twenty and a best fit of the logarithm of the numerator as a linear function of n, between 
n = 10 and n — 20, gives the value of A within an accuracy of 0.1%. For smaller values of A (typically A < 1.5) 
the number of iterations is larger, but the accuracy, for a given CPU time, is smaller. In such "difficult" cases one 
analytically finds the factorizations up to n — 7 and implement the first steps of these factorization schemes in the 
semi-numerical method. We are then almost guaranteed that no accidental factorizations will occur for n > 7, and 
therefore we can average over many initial matrices. Even so it remains difficult to discriminate between a truly 
polynomial growth jj],E| (A = 1) and an exponential growth with A ~ 1. The complexity values close to one clearly 
need to be revisited by other methods we present below. 

B. Equivalence relations between permutations 

Even if this semi- numerical algorithm is efficient it is quite time consuming to use it directly on the 9! permutations. 
To classify the complexities associated to a large set of (birational) transformations like the one associated to the 9! 
permutations of 3 x 3 matrices, one certainly needs to reduce this set as much as possible. For instance one can 
try to find symmetries such that two permutations, related by the symmetry, yield the same complexity A. These 
symmetries allow to build equivalence classes and, thus, to restrict the exhaustive analysis to a only one representent 
in each class. Furthermore one may have the prejudice that any non trivial symmetry could enable to explain a 
possible integrability structure of the mappings and beyond, the structures associated with the classification of the 
Arnold complexity of these mappings. 

There actually exist quite trivial symmetries, corresponding to relabeling of rows and columns JT4j , for which the 
complexities of the associated K 's are obviously equal. It is possible to go a step further and define a set of equivalence 
relations TZ^ between the permutations, yielding new equivalence classes such that any two permutations in the same 
"new" equivalence class, 1Z^ n \ automatically have the same complexity A. Equivalence relation TZ^ amounts to saying 



5 However one should keep in mind that there is nothing specific with 3x3 matrices. These results simply generalize to q x q 
matrices (see for instance M). 
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that two equivalent permutations are such that the n-th power of their associated transformations K are conjugated 
(via particular permutations, product of row permutations, column permutations and possibly the transposition, see 
appendix A for more details). An exhaustive inspection has shown that the equivalence relations IZ^'s "saturate" 
after n — 24: with obvious notations TZ^°°^ — TZ^ 24 ^ . One finds out that the "ultimate" TZS°°' equivalence classes 
can only have 72, or 144, elements. Among the "ultimate" 7^°°) classes, one wants to distinguish between the 
classes that were already classes, that we will denote from now on TZ^ , or IZ 1 ] 1 ^, according to their number of 
elements , and the other ones we denote 7£ 7 2 or 7Z^ . Being an 7^°°) equivalence class which does not reduce 
to a VS 1 ^ equivalence class, me ans the existence of several non trivial relations between the permutations in the 
7^(°°) equivalence class (see (A3) in appendix A). This implies strong constraints on the respective orbits. One thus 
expects more properties, and structures, inherited from this fact. The 362880 permutations are grouped into 2880 
equivalence classes (instead of 30462 "relabeling" equivalence classes in [[l4]). In Tab. | the number of the respective 
72-72 > ^-144) - ^-72°' ) or 7^144 classes is displayed. Since the complexities do not depend on the chosen representent, we 
picked a representent in each IZ^ 00 ^ class and performed, for it, the semi-numerical method previously explained. 

For 3x3 matrices, the complexities are necessarily such that : 2 > A > 1 . Remarkably, instead of getting a quite 
complicated distribution, or spectrum, of values for the complexities, we have obtained values which are always very 
close, up to the accuracy of the method, to a set of seventeen values given in the left column of Tab. | (see below) and, 
of course, the integrable value A = 1. To test the accuracy of the method we got complexities for two representants 
of the same class (that should, as we know, have the same complexity value exactly). We always obtained an equality 
of the corresponding complexities, up to an error of 10~ 3 . This accuracy is however not always sufficient enough to 
discriminate between some complexities displayed in the left column of Tab. |. In order to fix our mind it is necessary 
to obtain the exact expressions of these complexity values, for instance by getting the factorization scheme (^), (Q), 
and thus the generating functions a(x) and (3(x). 



C. Revisiting the complexity spectrum via exact factorization schemes 

For most of the TZ^ 00 ^ equivalence classes (2832 out of 2880), the complexity values, obtained with our semi- 
numerical method, are extremely close to the upper limit A = 2 . In fact one can figure out that these complexity 
values are actually exactly equal to 2. Therefore we can focus on the analysis of the remaining 48 classes, finding 
systematically their factorization schemes and associated generating functions. We actually found these factorization 
schemes and the associated generating functions, and were actually able to see that the previous numerical spectrum 
exactly corresponds to eighteen algebraic values listed in Tab. |[ Among these eighteen algebraic values, let us take four 
illustrative examples. We give for each example, the permutation representing the TZ^°°^ equivalence class, the value 



of A , and \i, defined in section IB, the expressions of (3(x) and p(x), since they respectively correspond to the simplest 
"degree generating function" and "exponent generating function" . The other generating functions can be deduced 
from these two, using linear functional relations (^), (J|) and (JsJ) between the generating functions [^). Furthermore 
relation ( |lO| ) remains valid for all the factorization schemes associated with all the various permutations studied here. 
We first give the permutation t itself, using the notation, already used in [Q, where PoPiP2PzPiPbPaP7Ps, means that 
(tM)i — M Pi , the entries of the matrix being enumerated consecutively (i.e. M\\ — Mo, M\<2 — Mi, M13 = M2, 
M 2 i = M 3 , M 33 = Ms)- 

• First example. Permutation 407326518 yields A ~ 1.61803 • ■ • and /i = 1 and : 

/?(*) = I"* 2 , x) = 1 (w) 

3 x 1 — x — x 2 ' (1 — x) 2 ■ (1 + x) 

• Second example. Permutation 417063582 yields A ~ 1.83928 • • • and fx ~ 1.32471 • • • and : 

0(x) _ 1 - x 2 - X 3 _ (1 + x) ■ (1 - x + z 4 ) 

3x (l-x) 2 -(l + x)-(l ~x - x 2 - x 3 ) ' P[X > ~ l~x 2 ~x 3 [ ' 

• Third example. Permutation 164273085 yields A ~ 1.83928 • • • and [i — 1 and : 

P(x) 1 + x + x 2 . , 1 + x 3 + x 4 + x 5 M „. 

~3x~ = l-x-x 2 -x 3 ' P{X) = 13^5 ( 16 ) 

• Fourth example. Permutation 174528603 yields A ~ 1.97458- ■ • and fx ~ 1.32471 ■ ■ • and : 
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(3(x) 1 - x 2 — x 3 1 - x + x 7 + x 8 

~Jx~ = (l-x)-{l-x-2x 2 -x 3 +x 4 + 2x 5 + x e ) ' P ^ X > = {1-x + x 2 ) ■ (l-x 2 -x 3 ) ( ' 

The exhaustive analysis of the factorization schemes, and the associated degree, and exponent, generating functions 
(a(x), (3(x), T](x), <f>(x) and p(x)), confirms that the complexities are actually independent of the representent in the 
equivalence class. On the contrary, the factorization schemes and the associated degree, and exponent, generating 
functions may depend^ on the chosen representent in the equivalence class. In other words to two permutations in the 
same class of equivalence correspond the same (up to 1 — x or 1 + x , or iV-th root of unity factors) denominators for 
the degree generating functions a(x), (3{x). By contrast the numerators, as well as the exponents generating functions 
are representent dependent (see the previous four examples). Most of time the stability of the factorization scheme 
and thus, in a second step, the occurrence of rational upnsraXAng functions, corresponds to a simple periodicity of the 
exponents rj n , 4> n or p n in the factorization scheme ([|), (ji|). This periodicity is simply associated to the fact that the 
exponent generating functions have N-th root of unity poles : 1 — x 2 , 1 — a; 8 , l — x e , ■■■ (see p(x) in ([To])). However 
one sees, on examples (|l5| ) and (17), that one may have a stability of the factorization scheme with an exponential 



growth of these exponents r\ n and <p n . These exponent generating functions, of course, have a "scheme complexity" 
p smaller that the growth complexity A. This "scheme complexity" p is the inverse of the poles of p{x), <fi(x) or 
n(x), that is (for (^)), p ~ 1.32471 • • • < A ~ 1.83928 • • ■. Recalling @, for which p = 1 and A ~ 1.8392 • ■ ■, and 
(|l5|), one sees that the same growth complexity A can be associated to several "scheme-complexity" p. Conversely, 
comparing the fourth example (17) and the second example (|l5|), one sees that one "scheme-complexity" p can actually 



yield several growth complexities A. 



D. To sum up 



All these factorization scheme calculations confirm the results of the semi-numerical method and are summarized in 
Tab. H Most of the 362880 birational transformations considered here do correspond to the most "chaotic complexity" , 
namely the upper bound A = 2 : one has 359568 such A = 2 birational transformations, that is 99.0873 % of 
all the birational transformations. It is known 14 1, that some symmetry-classes correspond to situations where the 
determinantal variables x n 's, defined by @, are periodic (denoted "Period." in the Table ||). This x n = x n+ N 
situation may correspond to situations where mapping K, itself, is of finite order (trivial integrability) , but also 
to polynomial growth situations, that is, A = 1 exactly. One remarks that IZ^^ contains all the integrable, or 
polynomial growth, mappings and, up to one class in 72^ , all the mappings such that x n — x n +jsr, including the 
situations where mapping K, itself, is of finite order. 



III. VARIOUS GENERALIZATIONS 



We now show that all these results also apply for a much larger set of rational transformations. The number of 
permutations of entries of 3 x 3 matrices being finite it has been possible to perform an exhaustive analysis. For more 
general transformations, depending on continuous parameters, is not anymore possible and we will proceed just with 
chosen examples. These examples always combine homogeneous transformations of the entries of a matrix together 
with the matrix inversion. Therefore the transmutation relations detailed in appendix A still apply, yielding again 
non-trivial symmetries for these new set of transformations. 



A. Combining different K's 



Let us first consider permutation 146237058 , and its associated A ~ 1.97481 • • • transformation K%, and permu- 
tation 471562380 and its A ~ 1.54258- •• transformation Ki- Let us compose the two previous transformations. 



Considering one TV 00 ' equivalence class, one does not get as many factorization schemes as the number of elements in the 
equivalence class. It seems, inspecting directly all the 9! factorization schemes (but only up to twelve iteration steps), that, 
most of the time, one gets, at most, two possible factorization schemes for a given 7?.'°°' equivalence class, and that the set 
of all the possible factorization schemes would be twenty one (besides the polynomial growth situations which can be quite 
"rich"). 



6 



A 


Polynomial 


'v 144 




/v 144 


;v 72 


Total 


Total 




2146 


660 


14 


60 


2880 


2 


1 - 2x 


2145 


640 
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33 


2832 
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1 - 2 x + x A - 2 x' 1 + x 4 - 2 x a + x° 





2 








2 


1.974584654 


1 - x - 2 x' - x J + X* + 2 x a + x u 





1 








1 
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1 -2x + x° -x' 





2 
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1.946856268 
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1 
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1 








1 
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1 - 2 x + x A - 2 x A + 2 x* - 2 x 3 











1 


1 


1.88320350 


1 - 2 x + x^ - 2 x' 1 + x 4 





2 





6 


8 


1.866760399 


1 - 2x + x J - x 4 





1 








1 


1.860073051 


1 - x - x^ - x 4 - 2 ■ x a 





1 








1 


1.857127516 


1 - 2x + x^ - x J - x° - x' + x B - 2x M + x IU 





1 








1 


1.83928675 


1 - x - x^ - x A 





2 








2 


1.75487766 


l-2x + x^-x J 


1 











1 


1.61803399 


1 — x — x" 3 





3 








3 


1.57014731 


1 - x - x J - x 3 





1 








1 


1.542579599 


1 - x - x J - x' - x B 





1 








1 


1.46557123 


1 + X - X J 











2 


2 


1 ( Pol.gr.) 


1 - x , 1 - x" , ■ ■ ■ 











9 


9 


1 ( Period.) 







1 
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10 



TABLE I. 



From these two "atoms" we build the "molecule" K, = Ki ■ K\. Note that K = K\ ■ K2, obviously has the same 
complexity. 

This example is an interesting one since the complexity (obtained from the previous semi-numerical calculations) 
of the "molecule" K. = K2 ■ K\ is smaller than the product of the two complexities of K\ and K2 : A(/C) ~ 2.897 < 
1.9748 • 1.5426 ~ 3.0463. In general the combination of two complexities Ai and A2 gives a complexity for the 
"molecules" larger than the product Ai • A2 , often equal to the upper bound (here A uppcr = 4). 

The factorization scheme of K, is of the same type as the one described in jl^l , namely a 11 parity- dependent' 
factorization scheme. It is detailed in appendix (|b|) and yields a degree generating function f3(x) : 

P{x) _ 1+ 2x - x 2 -x 4 + x 6 

~Jx~ ~ 1 - 3 x 2 + x A - x 6 - 2 a; 8 ^ ' 

The complexity of the molecule K. does not identify with the complexity of K%, or the one of Ki- It is a true 
new algebraic number. This algebraic expression for the complexity of the molecule is in good agreement with the 
semi-numerical value obtained above. We have systematically studied such "molecules" for a choice of eighteen 
representants of the eighteen complexities of table | combined with themselves, and beyond, with other representants. 
If one barters the permutation t^ for another representent t 2 in the same TZ^°°^ class, transformation K2 being 
modified accordingly (K2 — > ), the new "molecule" IC^ = ■ K\ yields, in general, another algebraic 

value for the complexity A : the equivalence relation IZ^ 00 ^ is no longer compatible with the "molecular structure". 

For all these "molecules" the parity- dependent factorization scheme, yields algebraic numbers for the complexities 
of these molecules in agreement with the values obtained from our semi-numerical now applied for the "molecules" . 
Combining among themselves all the permutations yields a large number of different algebraic complexities, much 
larger than the number of complexities obtained combining only representants of the 7^°°) classes among themselves. 



B. From permutations to linear transformations 

We got algebraic results on birational mappings associated with permutation of the entries. We now address 
the following question : are these structures (existence of a stable factorization scheme) dependent of the fact that 
we are dealing with permutations ? In other terms, does one loose these algebraic properties when deforming the 
permutations in most general transformations ? The most simple, and natural, generalization amounts to replacing 
the permutation of the entries by linear combination on the entries. 

Let us now consider a first example, namely the quite general linear transformation depending on twenty one 
parameters : 
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m M mi,2 mi, 3 

ni2,i ™2,2 rn 2 ,3 — ► (19) 
m 3) i to 3 , 2 m 3i3 . 

mi,i an mi,i + a i2 mi,2 + ai3 ^1,3 + «2i »^2,i + 022 "12,2 + 023 ™2,3 + a 3 i 7713,1 + a 32 m 3i2 + a 33 to 3j3 toi j3 

^2,1 C21 7772,1 + C22 m 2 .2 + C23 ^2,3 m-2,3 

m3,i 6n mi,i + 612 mi,2 + 013 mi,3 + 621 m 2 ,i + o 22 7772,2 + ° 2 3 ^2,3 + hi "13,1 + 032 m 3)2 + 6 33 m 3 , 3 ra 3 , 3 

This particular form singles out the rows of the 3x3 matrix (and thus can be understood as an RCT-compatible form, 
see appendix A). Similarly to the previous paragraphs let us introduce the homogeneous transformation K = L ■ I . 
Factorizations again occur at each iteration step. These factorizations correspond to a stable factorization scheme 
giving a growth like X N , where A ~ 1.61803 • • •. It is of the general type described in (||) and (Q). This yields the 
following generating functions : 



3 X 



1 



p(x) 



1 



(20) 



These results are actually valid for any "sufficiently generic" choice of the twenty one parameters. One thus has a first 
"universality" property: the complexity A is "generically" not dependent of the previous twenty one parameters. Fur- 
thermore relation (^|) (and consequently relation (|l5|)) remains also valid for all the factorization schemes associated 
with all the linear transformations studied in this section. Complexity A ~ 1.61803 • • • (corresponding to polynomial 
1 — x — x 2 ) is a complexity value already found in table |[ in the sixteenth row. It is noteworthy that no choice of the 
twenty one parameters leads to a permutation of any of the three classes corresponding to A ~ 1.61803 • • •. Besides 
the identity, the only choice of parameters, leading to a permutation, is &12 = a 3 2 = 1, all others being zero. The 
permutation is then the transposition^ Mi,2 <-> Af 3i 2 which corresponds to the complexity A ~ 1.46557- ■ •. This 
transposition is not isolated in the twenty-one parameters set of transformations. Actually if the parameters verify 
the two conditions 612 • a 3 2 — ax% ■ o 3 2 = 1 and ai2 + o 3 2 = , all the other parameters being zero, one then gets 
additional factorizations, the modified factorization scheme yielding : 



p(x) = 



1 



l-x' 



3x 



1 



1 -x 



(21) 



and again A ~ 1.46557- ■ -. There is also polynomial growth subcases, for instance 012 = C22 = 032 = 1, o-u and 022 
arbitrary non zero, all the other ones being zero. Generally speaking, having a complexity generically independent 
of r parameters (here twenty one) , one can only expect more factorizations on some subvariety of the r-dimensional 
space, and consequently a smaller complexity value A on this very subvariety. 

We now give another eleven parameter example associated with the following linear transformation : 



L 



7711,1 mi, 2 TOi )3 

m 2 ,i m 2l 2 m 2 ,3 
m 3 ,i m 3i 2 m 3>3 



"11,1 TOi, 2 + &21 TO 2 ,1 + 2 2 "7-2,2 + »23 W 2 , 3 ?™13 

777.2,1 012 mi,2 + a 2 i m 2 ,i + a 2 2 m 2 _ 2 + a 23 to 2 , 3 + a 32 m 3j2 m 2j3 

»T73,1 m 32 + C 2 1 7772,1 + C22 m 2:2 + C 2 3 7772,3 ^3,3 



(22) 



For K = L ■ I the corresponding generating functions are : 

1 0(x) 



p(x) 



1 



3 x 



1 



(23) 



The numerator of (3{x) does not appear in table Tab. |: this mapping has a new value for the complexity A ~ 
1.72088 • • ■, not previously obtained for any of the 9! permutations. 

Family (^2|), depending on eleven continuous parameters, also enables to address the following problem : is the 
complexity growth crucially dependent on the reversible character |l9| ] of the transformations ? In fact one may lose the 
birational character of K when, for instance, the linear transformation L becomes singular. This is very easy to realize 
for some condition on the eleven parameters (codimension one subvariety). For instance, taking 622 = 2 , 022 = 87 , 



7 Denoted class IV in [HE 
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a\2 = 5, a 3 2 = 7, C22 = 11 , all the other parameters being zero, leads to a non invertible mapping K = L ■ I. 
One easily verifies that the factorization scheme, the associated generating functions and thus the complexity A, 
are unchanged in this case and, more generally, on such singular subvarieties. With this first rational, non invertible, 
example one sees that the rational character of the generating functions is not a consequence of a "simple" invertibility 
of the mapping (see also Q ) . 



C. From linear transformations to homogeneous polynomial transformations 



There is nothing specific with linear transformations. For instance, let us consider the following quadratic trans- 
formation depending on twenty one parameters (which is reminiscent of (|l9|)) : 



mi,i mi, 2 mi, 3 

7712,1 1712,2 "T-2,3 

m 3 ,i m 3 . 2 m 3:3 

m\ x an m\ x + a 12 m\ 2 + a 13 mf 3 + a 2 i m\ x + a 22 m\ 2 + a 23 m\ 3 + a 3i m\ x + a 32 m| 2 + 033 rn% 5 m\ 3 



(24) 



c 2 i m\ x + c 22 ml 2 + C23 rn\ 3 



'•2,3 



n\ x bn m\ Y + bi2 m\ 2 + 6 13 m 2 3 + b 2 \ m\ Y + 622 m 22 + 6 23 m\ z + 6 3i m% x + 6 32 m\ 2 + b 33 m 3 3 mf g 

The homogeneous transformation A' = Q ■ I gives again a stable factorization scheme. In this case, where Q 
is no longer a linear transformation, but a homogeneous polynomial transformation of degree r (here r = 2), the 
factorization scheme remains of the general form (^) and (|]). As far as generating functions are concerned some 
modifications have to be done. Firstly the p n 's, and associated p(x), should be replaced by the 7„'s defined by : 



K(M n ) 



K{M n ) 
det(M n y 



M, 



n+l 



f7o f7l fl 
Jn+1 Jn ' Jn 



(25) 



and t he c or resp onding generating function j(x). The linear relations between rj(x), 4>(x) and j(x) are slightly modified 
(see (|ci|), ( |C4[ ) in appendix C). Secondly, relation ( p^0| ) is no longer valid here. A new relation has to be introduced 
playing the same role. Transformation K — Q ■ I is a homogeneous transformation of degree — r. Instead of 
introducing the determinantal variables x n through (^J), let us introduce x n by : 



x n (M ) = det(K n+1 (M ))-(det(K n (M ))Y 

These new determinantal variables x n are well-suited ones since they are invariant under a rescaling of Mq 
x n (Cst ■ Mq) = x n (M ) ■ Relation ( |l0| ) becomes : 



Wi fW 2 



rW 3 



f 

Jo 



(26) 



(27) 



Again one can introduce the generating function of these exponents W n and see that relation (13) still holds. From 
the stable factorization scheme of K = Q ■ I one now gets : 



0(x) 



1 



3x l-3z-2x 2 ' 7(:E) 2 'l-2x 

This gives a complexity value A ~ 3.5615 • ■ • . Let us consider the expression of a(x) 

3 • (1 +x - 2x 2 + 4x 3 ) 



1 



a(x) — 



(1 + 2x) ■ (1 - 2x) • (1 - 3x - 2x 2 ) 



(28) 



(29) 



On this expression one sees that other poles occur. The inverse of these additional poles, namely ±2, are actually 
smaller that the complexity value 3.56155 • • • . The existence of "subdominant" poles already occurred with permu- 
tations of entries, or linear transformations (see /3(x) in (|l5|)): we often had 1 — x, or 1 + x, additional factors in 
the expressions of the degree generating functions. With expression (^) one sees the occurrence of a 1 — 2 -x factor 
instead of 1 — x factors. 

There is also nothing specific with quadratic transformations. Let us introduce the simple homogeneous polynomial 
of degree r : 
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mi, 2 


™1,3 " 




' m l,l 


rrf l2 


m l,3 " 


Qr ■ 


m 2 ,i 


™2,2 


™2,3 




m 3,2 


m 2,2 


m 3,l 




. TO3,1 


™3,2 


™3,3 - 




_ ™2,3 


m 2,l 


m 3,3 . 



and its associated homogeneous transformation K 
r > 2 (for r = 1 transformation (^0|), and = Q r 

K{M n -i) 



(30) 



= Q r ■ I . Its factorization scheme is very simple, it reads for 
/, become trivial) : 



M n = 



fr 

Jn-1 



det(M n ) = /„ +1 • f n 



which yields the following linear relations on the a„'s and /3„'s (see also appendix C) : 

a n = 2 • r • a„_! - 3 • r ■ /3 n _i , a n = (} n+1 + 2-/3 n 

It gives the following generating functions for arbitrary r > 2 : 

0(x) _ 1 



3a; 



1 + 2- (1-r) -x 



rj{x) — r, 4>( x ) = l + 2>x. 



7(x) 



(31) 



(32) 



(33) 



For homogeneous polynomials of degree r one can show that subdominant poles, like 1 —r-x, may occur instead of 
the previous 1 — x and 1 — 2 x factors. 

For r — 2 , one remarks that one gets a degree generating function : 



3 x 



1 



1 - 2 a; - 2 a; 2 



(34) 



which is not the limit of (|28|). The generic complexity corresponding to (|2q), namely A ~ 3.56155 • ■ • is changed, 
for (|3(]) taken for r — 2, into A ~ 2.73205 • ■ -. There actually exist many subvarieties of the twenty-one parameter 
space of transformations ( |24| ) on which the generic complexity A ~ 3.56155 • • ■ is modified into another (smaller) 
algebraic value. One remarks that the subvarieties of the twenty-one parameter space of transformation ( |l9| ) (for 
instance, at2 — C22 = &32 = 1, an and 022 arbitrary non zero, all the other ones being zero, previously mentioned as 
a polynomial growth subcase) also yield "non-generic" complexities for (p3). 



IV. CONCLUSION 

In previous papers Jl7],^3| it has been shown that the topological entropy, and the Arnold complexity, actually iden- 
tify on various simple two-dimensional birational examples, and that these quantities are actually algebraic numbers. 
The generating functions corresponding to these two complexity measures, namely the dynamical zeta function l2p|-|2^1 
and the various "degree" generating functions (like /3(x)) were shown to be simple rational expressions with integer 
coefficients jf?]], the dominant poles in these two sets of generating functions being the same. When one analyzes 
birational transformations depending on more than two variables, it becomes very difficult to calculate even the first 
coefficients of the expansion of the dynamical zeta function. On the contrary the calculations on the degree generating 
functions (associated to the Arnold complexity) can be quite easily performed, even for birational transformations of 
many variables (the q 2 entries of a matrix 0). 

Analyzing exhaustively a first finite set of 362880 birational transformations (associated with all the permutations 
of 3 x 3 matrices), we have obtained non-trivial, but still simple, "spectrum" of eighteen algebraic Arnold complexities 
for the corresponding dynamical systems. In a second step it has been shown that these results can be drastically 
generalized along three different lines preserving the algebraic character of the complexities. Firstly, one can com- 
bine these birational transformations together, and get extremely rich sets of algebraic complexities. Secondly, one 
can consider (generically birational) transformations, associated with linear transformations of the entries of 3 x 3 
matrices, and still get sets of algebraic complexities. Remarkably one has another universality property here: these 
algebraic complexities do not depend on many of the continuous parameters associated with the linear transforma- 
tions. Thirdly, one still gets sets of algebraic complexities with rational transformations (associated with homogeneous 
polynomial transformations on the entries) which again can depend on many continuous parameters. With this last 
generalization we have completely lost any invertible character of the transformations. On the top of that these 3x3 
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matrix calculations can be simply generalized to q x q matrices for arbitrary^ q. Combining several of these rational 
transformations depending on several continuous parameters together, one certainly gets again rich sets of algebraic 
complexities. 

We end up with an extremely large set of transformations, so large that, clearly, it should be a powerful tool to 
study discrete dynamical systems. 

APPENDIX A: A "TRANSMUTATION" PROPERTY OF THE MATRIX INVERSION 

Let us sketch here some non trivial symmetries between the permutations. 



The transformations, considered in sections ( |IB| ), (II C), are products of matrix inversion and permutations of 



the entries. Any such non trivial symmetry of the birational transformations K should correspond to a non-trivial 
relation between matrix inversion and permutations of the entries of the matrix. Such relations actually exist. They 
correspond to a "transmutation" property between the inversion and permutations P and Q. There actually exist 
two permutations P and Q such that : 

pi = Tq (ai) 



Permutations, such that a "transmutation" relation (AI) is satisfied, do exist : one can easily build examples 



by combining product of permutations that permutes only rows of a q x q matrix (that we will denote by "-R"), 
permutations that permutes only columns of a q x q matrix (that we will den ote by "C") and, possibly, the matrix 
transposition we denote "T". Examples of permutations P and Q, such that ( |Al[ ) is satisfied, read : 

P = R-C -T e , where : e = , or 1 (A2) 



and similarly for permutation Q. A permutation P having such a decomposition (A2) will be called an "RCT" 
permutation. 

Let us consider two permutations t\ and t2, yielding respectively the two birational transformations K\ = t\ ■ I 
and K2 — t<i ■ I ■ Let us introduce the following relation of equivalence between two permutations t\ and ti ■ t\ and 
t2 will be related if they are such that there exists an "RCT" permutation, bo, such that : 

K? = bv-KZ-bv 1 (A3) 

Relation (|a|) can easily be seen to define a relation of equivalence between t\ and we will denote IZ^ : 

h Tl {n) t 2 (A4) 

Note that this VS n ^ equivalence relation is compatible with the inverse in the permutation group t — > t^ 1 . 
Also note that the equivalence of two permutations, up to simple rows and columns relabeling, is an equiva- 
lence, however, conversely, the 7Z^ equivalence does not reduce to the simple, and quite trivial, equivalence of two 
permutations up to simple rows and columns relabeling. Obviously rows and columns relabeling of the matrices do 
not modify their integrability properties fL4| , as well as the growth of the calculations. 

It is obvious that if t\ TZ^ t% then t\ 7?/ nxp ) t 2 f° r an y natural integer p. This is a consequence of the fact 
that : 

K[ l = b Q ■ K% ■ b- 1 yields : K™ p = b Q ■ K^ p ■ b' 1 (A5) 

If two permutations, t\ and are in the same equivalence class with respect to VS m \ and if £2 and £3 are in the 
same equivalence class with respect to TZ^ where n ^ m, t\ and £3 are in the same equivalence class with respect 
to 7^(™ xm ) j or with respect to VS- N ^ for some "large enough" integer N. In fact it can be shown, on the example 
of the equivalence classification of the permutations of 3 x 3 matrices, that this value of N corresponding to the 
( "asymptotic" equivalence) relation is actually equal to N = 24 . 

If two permutations t\ and t<i are in the same equivalence class, with respect to TZ^ m \ the complexities (which are 
real positive numbers), associated wi th their respective birational transformations K\ and K 2, we denote Ai and A2 



are, as a straight consequence of (A3), related by 



8 The "spectrum" of values of the complexity A depends on q, see for instance Q. 
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A™ = A™ (A6) 

Therefore one sees that their complexities are equal : Ai = A2. In particular if one considers the (largest) equivalence 
classes corresponding, for 3x3 matrices, to TZ^ 24 \ all the representants in one of these V} 24 ^ equivalence classes will 
have the same complexity growth A. 

APPENDIX B: A "MOLECULAR" FACTORIZATION SCHEME 

The factorization scheme of K. = t\ ■ I-t-z -I, corresponding to permutation 146237058 and permutation 471562380 



(see section (III A), is of the same type as the one described in |17|, namely a parity- dependent factorization scheme 



(which is a straight consequence of the fact that one acts with Ki, and then with K2, and again ...) : 

h = det(M ), Mi = Kt(M ), f 2 = det(Mi) , M 2 = K 2 (M{) , / 3 = det ^) ) Ma = Kl (M 2 ) 

J2 



t J*tM\ M w(iJt\ t detjMj) KM) det(M 5 ) 

/ 4 = det(M 3 ), M A = K 2 (M 3 ), f 5 = 3 , M 5 = , / 6 = , 2 , ■•• (Bl) 

J 2 ' U h J 2 ' U 



and for arbitrary n > 3 



det(M n ) — f n+ i ■ f n ■ fn- 2 ■ fn-6 ' fn-8 ' /n-10 ' /n-12 ' /n-14 ' 1 ' 

#i(M n ) = M n+1 -/ n _ 2 (B2) 



for n even and 



det(M n ) - f n+1 ■ / n _i • / n _ 3 • / n _ 5 • / 2 _ 7 • f n _ g ■ f n _ n ■ / n _ 13 • ■ ■ 

Jf 2 (M n ) = M n+1 • /„_ 3 • /„-7 • fn-9 ■ fn-U • /n-13 ' ' ' (B3) 

for n odd. This yields the following expressions for the odd and even parts of a(x) and (3(x) ("2" for even and "1" 
for odd) : 

Bo(x) - 6 - x2 bax) - 3-^a + ^ 2 )-(-i+^) 2 (^+i) 2 

P2W 1 -3a; 2 +.t 4 - a; 6 -2a; 8 ' Pl[ ) 1 -3x 2 +x 4 -x 6 -2x 8 

3-(l +4a; 4 -4a; 6 +a; 8 ) 6 ■ x ■ (1 + x 4 - x % + a; 8 ) 

a2 ^ X ' ~ (l-a; 2 )-(l -3a; 2 + x 4 - x e -2a; 8 ) ' ~ (1 - x 2 ) • (1 -3a; 2 + x 4 - x e -2x 8 ) ( ' 

These generating functions yield a "molecular complexity" : A ~ 2.8581 • • •. These generating functions verify a 
parity dependent system of functional relations which generalizes the ones described in @] : 

x-ai(x) - (3 2 {x) = F 2p (x) ■ f3 2 (x) , x ■ a 2 (x) - (3 x (x) = F lm (x) ■ f3 2 (x) , (B5) 

a 2 (x) - 3 — 2 ■ x • ai(x) + 3 • G 2p ■ P 2 (x) = 0, oti(x) - 2 ■ x ■ a 2 (x) + 3 • G\ m ■ (fa(x) = 

where : 

F 2p {x) = x 2 + 2x 4 + x 6 + ^— F lm = 2x 3 -x 5 + -^ T , G Xm {x) = x 3 , G 2p = x 4 + 

1 — ar 1 — x l 1 — x £ 



APPENDIX C: EXPONENT GENERATING FUNCTIONS FOR HOMOGENEOUS POLYNOMIAL 

TRANSFORMATIONS OF DEGREE R 

Let us consider a homogeneous transformation Q r of degree r (like (p0|), or like (|24|) for r = 2) and its associated 
homogeneous transformation K = Q r - 1 . Relations (^) , (Q) are still valid but yield a slight modification of the linear 
functional relations @ and ^ , namely : 

((q - 1) • r ■ x - 1) • a(x) + q - q ■ x ■ r)(x) ■ j3{x) = 0, (CI) 
x ■ a{x) = (j){x) ■ (3{x) (C2) 
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Let us recall that, for ho mogen eous transformations of degree r, one must introduce, instead of p(x), the generating 
function 7(3;) (see section ( [II C )) defined by : 



K{M n ) 



K{M n ) 



+ 1 



det(M n y f n ° +1 ■ Pn 1 ■ fnU ■ 
This last relation yields a new relation : 

q + q ■ j(x) ■ (3{x) = (1 + r ■ x) ■ a(x) 
which has to be compatible with the previous two flC]] ), ( |C^ ) : 

r ■ 4>{x) — 7(2;) + x ■ r)(x) 



(C3) 

(C4) 
(C5) 
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